Predator cannibalism can shift prey community composition toward dominance by small prey species

Abstract Cannibalism among predators is a key intraspecific interaction affecting their density and foraging behavior, eventually modifying the strength of predation on heterospecific prey. Interestingly, previous studies showed that cannibalism among predators can increase or reduce predation on heterospecific prey; however, we know less about the factors that lead to these outcomes. Using a simple pond community consisting of Hynobius retardatus salamander larvae and their associated prey, I report empirical evidence that cannibalism among predators can increase predation on large heterospecific prey but reduce that on small heterospecific prey. In a field‐enclosure experiment in which I manipulated the occurrence of salamander cannibalism, I found that salamander cannibalism increased predation on frog tadpoles but reduced that on aquatic insects simultaneously. The contrasting effects are most likely to be explained by prey body size. In the study system, frog tadpoles were too large for non‐cannibal salamanders to consume, while aquatic insects were within the non‐cannibals’ consumable prey size range. However, when cannibalism occurred, a few individuals that succeeded in cannibalizing reached large enough size to consume frog tadpoles. Consequently, although cannibalism among salamanders reduced their density, salamander cannibalism increased predation on large prey frog tadpoles. Meanwhile, salamander cannibalism reduced predation on small prey aquatic insects probably because of a density reduction of non‐cannibals primarily consuming aquatic insects. Body size is often correlated with various ecological traits, for instance, diet width, consumption, and excretion rates, and is thus considered a good indicator of species’ effects on ecosystem function. All this considered, cannibalism among predators could eventually affect ecosystem function by shifting the size composition of the prey community.


| INTRODUC TI ON
Predators can strongly influence the abundance of their heterospecific prey, eventually affecting species richness and composition of prey community and its ecosystem functions. Population size (i.e., density) and foraging behavior of predators are fundamental elements in determining the predatory effects on their heterospecific prey (Ohgushi, 2012;Werner & Peacor, 2003;Wootton, 1994).
Notably, because individuals constituting a focal predator population can differ in their foraging behavior (Bolnick et al., 2011;Miller & Rudolf, 2011;Violle et al., 2012), the net effect of the predator population on their prey community would be the sum of each individual's effects. Therefore, studying factors affecting predator density and individual foraging behavior and how the changes in predator density and behavior are translated into the predatory effects will provide knowledge to better understand the spatial and temporal variation in community composition (Ohgushi, 2012;Werner & Peacor, 2003;Wootton, 1994).
Particularly for predators, cannibalism is a key intraspecific interaction affecting population density but also differentiating foraging behaviors of individuals constituting the population (Claessen et al., 2004;Fox, 1975;Polis, 1981). However, predicting how cannibalism among predators modifies their effects on heterospecific prey is challenging because cannibalism-induced changes in density and foraging behavior of predators can have contrasting impacts on the strength of predatory effects. In general, cannibalism among predators causes the following changes in their population. Firstly, cannibalism among predators reduces their density (Fox, 1975). Secondly, when cannibalism occurs, potential victims (i.e., non-cannibals) reduce their foraging activity to reduce their likelihood of being detected by the cannibals (Rudolf, 2012;Sih, 1982;Wissinger et al., 2010). Moreover, especially in predators in aquatic ecosystems such as fish, amphibian larvae, and aquatic insects, cannibalism often causes enhancement of growth of individuals that succeeded in consuming conspecifics (i.e., cannibals; Hardie & Hutchings, 2014;Kishida, 2011;Sniegula et al., 2019). The reductions in both density and foraging activity of non-cannibals are likely to reduce predation on heterospecific prey. In contrast, the enhancement of the growth can result in the improvement of the foraging ability of cannibals, thereby increasing predation on heterospecific prey. Notably, reductions in non-cannibals' density and foraging activity, and enhancement of cannibals' growth occur together within a predator population due to cannibalism. Therefore, the net effects of predator cannibalism on the strength of predation on heterospecific should differ depending on the relative importance between reductions in non-cannibals' density and foraging activity and improvement of cannibals' foraging ability. For example, cannibalism among predators should reduce predation on their heterospecific prey if the reductions in non-cannibals' density and foraging activity are more significant than cannibals' improved foraging ability in determining the strength of predation. Indeed, previous studies examining how cannibalism among predators affects their heterospecific prey showed both reductions (Crumrine, 2010a;Persson et al., 2003;Rudolf, 2006Rudolf, , 2012 and intensifications  of the strength of predation on the prey as a result of predator cannibalism. The differences in the reported patterns among studies raise a question; what factors affect the relative importance between the reductions in non-cannibals' density and foraging activity and the improvement of cannibals' foraging ability in determining the strength of predation on heterospecific prey?
The body size of heterospecific prey can play a pivotal role in determining the relative importance between the reductions in non-cannibals' density and foraging activity and the improvement of cannibals' foraging ability for the following reasons. Firstly, predators often co-occur with different-sized prey species, and also the prey community size composition varies across space and time.
Secondly, particularly for carnivorous fish and amphibian larvae, which are often top predators in freshwater food webs and exhibit cannibalism (Fox, 1975), the upper limit of consumable prey size is often restricted by the predator's gape size (Montori et al., 2006;Nosaka et al., 2015). It is worth mentioning that the upper limit of consumable prey size increases with the growth of predator individuals (Werner & Gilliam, 1984). All this considered, the direction of the net effect of predator cannibalism on the strength of predation on heterospecific prey can differ depending on whether focal prey are larger or smaller than the upper limit of consumable prey size for non-cannibals.
For example, when focal prey are larger than or similar to the upper limit of consumable prey size for non-cannibals, non-cannibals rarely consume the prey. Thus, the strength of predation on the large prey is less likely to be affected by the reductions in non-cannibals' density and foraging activity. Meanwhile, if cannibals reach large enough size to consume the large prey easily, the strength of predation on the large prey should be affected by the growth enhancement of cannibals. If this is the case, the strength of predation on the large prey should be intensified by cannibalism among predators eventually. On the other hand, when focal prey are smaller than the upper limit of consumable prey size for non-cannibals, both cannibals and non-cannibals consume the prey. Thus, the strength of predation on the small prey is likely to be affected both by reductions in non-cannibals' density and foraging activity and also by the improvement of cannibals' foraging ability. Notably, only a few individuals within a population can cannibalize and grow into gigantic cannibals (Huss et al., 2010;Kishida, 2011). All this considered, reductions in non-cannibal's density and foraging activity are likely to be more significant than the improvement of cannibals' foraging ability in determining the strength of predation on the small prey. Therefore, the strength of predation on the small prey should be reduced by cannibalism among predators eventually. Moreover, when the focal prey community includes both small and large prey, cannibalism among predators can further reduce the strength of predation on the small prey if cannibals shift their diet from small to large prey. In support of the hypothesized importance of prey body size in determining the direction of predator cannibalism effects, previous studies showing that cannibalism among predators reduces predation on prey used predator-prey systems in which focal prey items were smaller | 3 of 11 TAKATSU than the upper limit of consumable prey size for non-cannibals (e.g., predatory fish and their prey zooplankton; Persson et al., 2003).
Also, previous studies showing that cannibalism among predators increases predation on prey used a system in which the focal prey item was larger than or similar to the upper limit of consumable prey size for non-cannibals (predatory salamander larvae and their prey frog tadpoles; . However, so far, no studies have directly tested the hypothesis that cannibalism among predators can reduce predation on small heterospecific prey but increase predation on large heterospecific prey. The simple pond community consisting of the Hynobius retardatus salamander larvae and their prey Rana pirica frog tadpoles and aquatic insects (e.g., mayfly and chironomid larvae and water boatman), which is commonly observed in small ponds (one to several tens m 2 ) in Hokkaido, Japan, is an excellent model system to examine how cannibalism among predators affects the size composition of prey communities, because of the following knowledge. Firstly, salamander larvae consume their prey items, including conspecific victims (i.e., cannibalism), only when their gape size is larger than the body size of the prey (Nosaka et al., 2015). Secondly, there are significant differences in body size between R. pirica frog tadpoles and aquatic insects ( Figure 1). In early spring, both H. retardatus salamanders and R. pirica frogs lay their eggs in the small ponds. Several weeks after the end of the reproductive season of the amphibians, the recruitment of larval aquatic insects occurs. Due to the differences in the reproductive timing as well as species-specific body size, frog tadpoles, and salamander larvae are generally similar in their body size (i.e., frog tadpoles are large prey for the salamander larvae), while aquatic insects are far smaller than the larval amphibians (i.e., aquatic insects are small prey for the salamander larvae; Figure 1). For this study, I hypothesized that cannibalism among salamander larvae increases predation on frog tadpoles (large prey) but reduces predation on aquatic insects (small prey) at the same time.
Here, I report results from a field-enclosure experiment in which I manipulated the size structure of salamander larvae to control the occurrence of their cannibalism and then assessed how the effects of salamander cannibalism differ depending on differently sized prey items. The intervals between enclosures were more than 50 cm ( Figure 2a).

| Experimental environments
An emergence trap was placed on the top side of each enclosure to prevent metamorphs of amphibians and emerging aquatic insects from escaping from the enclosures (Figure 2a). Periphyton grows on the mesh wall as resources for frog tadpoles and aquatic insects.

| Experimental setting
A schematic diagram of the field experiment is shown in Figure 2b.

| Prey community
A prey community consisting of frog tadpoles and aquatic insects was established using the following methods. For frog tadpoles, I assigned 80 2-week-old laboratory-reared frog tadpoles into each of the 36 enclosures on 23 June 2013 (i.e., day 1 of the experiment; Figure 2b). Mean ± SD body length (snout-vent length) and body width of the frog tadpoles were 12.4 ± 0.8 mm and 7.7 ± 0.5 mm, respectively (N = 20) (Figure 1b). The density of the frog tadpoles (222 individuals m 2 ) was within their natural range (Michimae, 2006).
Methods for collecting and maintaining the frog tadpoles are described in Appendix A.
Since newly hatched aquatic insects can freely pass through the mesh wall of the enclosure when they are smaller than the mesh size (i.e., <1 mm), I used those which naturally colonized the enclosures for the experiment. While I observed several aquatic insects in the enclosures during the experiment, more than 97% of observed aquatic insects over the experimental period were composed of Ephemeroptera, Diptera, and Hemiptera (Appendix B). Hereafter, I focused on the Ephemeroptera, Diptera, and Hemiptera as focal aquatic insect prey in the experiment. Mean ± SD body length (from the tip of the head to the end of the abdomen) of Ephemeroptera larvae, Diptera larvae, and Hemiptera were 5.9 ± 1.2 mm (N = 20), 10.1 ± 5.1 mm (N = 20), and 3.8 ± 1.7 mm (N = 5), respectively ( Figure 1b). Mean ± SD body width of Ephemeroptera larvae, Diptera larvae, and Hemiptera were 1.3 ± 0.2, 1.2 ± 0.7, and 1.8 ± 0.7 mm, respectively ( Figure 1b). Therefore, those focal aquatic insects were far smaller than frog tadpoles ( Figure 1b). These indicated that aquatic insects were small prey and the frog tadpoles were large prey for the salamander larvae. In addition to the aquatic insects, I also observed zooplankton, such as copepods, in the experimental pond. The recruitment of such zooplankton provided additional prey for the salamander larvae during the experimental period. However, I did not include zooplankton as focal prey in this study because zooplankton were smaller than the mesh size of the enclosures (1 mm) during their lifetime and were too difficult to count.

| Experimental treatments
To control cannibalism among salamander larvae, I followed the methods used in the previous studies Takatsu et al., 2017). Since the occurrence of cannibalism among salamander larvae depend greatly on size asymmetry between interacting individuals , I manipulated size differences among salamander larvae (i.e., presence and absence of early-and latehatched salamander larvae) while keeping the total initial density of the salamander larvae constant across the treatments (i.e., 30 individuals, see below). In the experimental setting, the early-hatched salamander larvae can grow until the late-hatched salamander larvae hatch. This resulted in size asymmetry between them and eventually cannibalism (i.e., early-hatched salamander larvae can cannibalize late-hatched salamander larvae). I obtained the early-and late-hatched salamander larvae by manually controlling the water temperature experienced by the embryos from a single egg cluster (i.e., half-sib). Methodological details to obtain early-and late-hatched salamander larvae are the same as those shown in Appendix A in Takatsu and Kishida (2015).
Using the 36 enclosures, I established the following four treatments: (1) Cannibalism, (2)  was 6 days. The salamander larvae were assigned to the relevant treatments 1 day after they hatched (i.e., days 1 and 7; Figure 2b).
Mean ± SD body length and gape width of the salamander hatchlings at the assignment timing were 11.2 ± 0.8 and 2.8 ± 0.3 mm (N = 20), respectively. Methods for collecting and maintaining the salamander larvae are described in Appendix A. To make a No-salamander control treatment, I did not assign any salamander larvae to the remaining six enclosures. I adopted the unbalanced replication design to avoid excessive use of the animals. This is because previous studies using the same larval amphibian system showed that variances of the demographic and trait level consequences were larger in the Cannibalism treatment than in the No-cannibalism treatment Takatsu et al., 2017). Moreover, previous studies using the system showed that the mortality of the frog tadpoles in the absence of salamander larvae was very low (Kishida et al., 2014;Takatsu & Kishida, 2020). Each replicate was randomly assigned to one of the 36 field enclosures. The density of the salamander larvae (83 individuals m 2 ) was within their natural range (Michimae, 2006).

| Statistical analysis
All statistical analyses described below were conducted using R (ver-  (Table B2 in Appendix B).

| RE SULTS
There were no significant effects of salamander hatch timing difference (No-cannibalism -early vs. -late) on any of the focal demographic parameters (see Appendix C). Hence, the data from the two No-cannibalism treatments were pooled before conducting the following statistical analyses (hereafter, No-cannibalism treatment).  Figure  D2 in Appendix D]). Body length of individuals with the largest body length in the Cannibalism treatment (33.1 ± 3.2 mm) was 1.7 times larger than that in the No-cannibalism treatment (19.1 ± 1.3 mm) (ttest, t 28 = 17.35, p < .0001) (Figure 1b). Gape width of the largest salamander in the Cannibalism treatment (9.7 ± 1.3 mm) was 1.7 times larger than that in the No-cannibalism treatment (5.7 ± 0.6 mm) (t 28 = 11.70, p < .0001) (Figure 1b).

F I G U R E 3
Survivorship of the frog tadpoles during the experimental period. The thick horizontal bars represent the median, the box contains 50% of the data, and the whiskers indicate the range. Treatments not sharing the same lowercase letter were significantly different from each other (Tukey's HSD: p < .05). See Figure 2 Table 1). There was no significant effect of an interaction between treatment and census timing (p > .91) except for Ephemeroptera ( 2 2 = 50.61, p < .0001) ( Table 1). Thus, except for Ephemeroptera, I performed GLMM without considering the interaction effects, and  (Figure 4d,e). Posthoc pairwise comparisons found differences among treatments which are consistent with the patterns observed in the former analyses. Total number of emerging Ephemeroptera in the Nocannibalism treatment (5.2 ± 3.0 [mean ± SD]) was 38% and 84% lower than the Cannibalism (8.4 ± 4.0; z-ratio = 2.37, p = .047) and No-salamander (32.3 ± 9.4; z-ratio = −9.94, p < .0001) treatments, respectively. Total number of Ephemeroptera metamorphs in the Cannibalism treatment was also 74% lower than the No-salamander treatment (z-ratio = −8.24, p < .0001). Similarly, total number of emerging Diptera in the No-cannibalism treatment (352.2 ± 74.9 [mean ± SD]) was 17% and 28% lower than the Cannibalism (424.6 ± 84.8; z-ratio = 2.58, p = .027) and No-salamander (486.0 ± 64.6; z-ratio = −3.62, p = .0009) treatments, respectively.
Total number of emerging Diptera in the Cannibalism treatment was 13% lower than the No-salamander treatment while the difference between treatments was not significant (z-ratio = −1.53, p = .28).
The differences among treatments consistently observed in aquatic insects remaining in the enclosures and total number of emerging insects suggest that cannibalism among salamanders reduced predation on aquatic insects.

| DISCUSS ION
While cannibalism among predators can both reduce and increase predation on heterospecific prey (Persson et al., 2003;Rudolf, 2006Rudolf, , 2012, we know less about the factors that lead to these contrasting outcomes. Here, I tested the hypothesis that cannibalism among predators can reduce predation on small heterospecific prey but increase predation on large heterospecific prey using H. retardatus salamander larvae and their associated prey community. In the focal prey community, aquatic insects were small prey for salamander larvae, and frog tadpoles were large prey (Figure 1).
Consistent with my hypothesis, in the field-enclosure experiment, while the presence of salamanders reduced the abundance of the small and large prey, the numbers of the aquatic insects remaining in TA B L E 1 Summary of results of analyses examining the effects of treatment, census timing, and interaction between them on aquatic insects remaining in the enclosures using generalized linear mixed models (GLMM) with quasi-Poisson error distribution (log-link) Enclosure ID was included as a random factor to give a repeated measured design.
the enclosures and emerging aquatic insects were generally lowest when salamanders did not cannibalize, but survivorship of the frog tadpoles was lowest when salamander cannibalized (Figures 3 and   4). This study provides empirical evidence that the effects of predator cannibalism on prey depend on the prey body size.
There is, however, a possibility that the differences in the direction of predator cannibalism effects among prey items could be solely explained by prey phenology rather than prey body size. In There are two possible mechanisms causing the prey body size dependency in the effects of salamander cannibalism. Firstly, aquatic insects were within the consumable size range of both noncannibals and cannibals, although frog tadpoles were too large for non-cannibals to consume (Figure 1b). As a result, the reductions in non-cannibals' density and foraging activity due to cannibalism were less likely to affect the predation of the frog tadpoles. Therefore, the strength of the predation on the frog tadpoles was solely determined by the effects of the improvement of cannibals' foraging ability, as shown in the previous studies . At the same time, the predation on the aquatic insects was likely to be affected by reductions in non-cannibals' density and foraging activity and improvement of cannibals' foraging ability. However, possibly because only a few individuals within a population can cannibalize and become giant cannibals (Kishida, 2011), the effects of improvement of the foraging ability of the few cannibals cannot overwhelm the effects of reductions in non-cannibals' density and foraging activity. Indeed, there were less than five salamanders whose gape size was larger than the body width of frog tadpoles, while cannibalism reduced the population by less than half ( Figure Figure G1 in Appendix G). Therefore, it is expected that the relative size relationships between salamander larvae and their prey were more important than the diet broadening of cannibals in explaining the contrasting effects of salamander cannibalism on predation.
Considering that predators often broaden their diet and or shift their diet from small to large prey with their growth (Werner & Gilliam, 1984), paying attention to the size composition of prey communities is likely to be essential when predicting communitylevel consequences of predator cannibalism. For example, suppose that a focal prey community is dominated by large prey species (e.g., frog tadpole in this study system). Then, cannibalism among predators can shift the prey community toward dominance by small prey species by increasing predation on large prey but reducing predation on small prey simultaneously. Importantly, as well as species identity, the body size is considered as an important factor in explaining variation in individuals' ecological roles since body size is correlated with, for instance, trophic position, consumption rate, and nutrient excretion rate (Romanuk et al., 2011;Vanni et al., 2002;Werner & Gilliam, 1984;Woodward et al., 2005). Partly due to this, large species often play a disproportionately large role in determining ecosystem functions, such as primary production (Séguin et al., 2014). All this considered, the shift in the prey community toward dominance by small prey species due to predator cannibalism could have significant impacts on ecosystem function. To further our understanding of community-level consequences of predator cannibalism, examining how predator cannibalism effects change along with the size composition of prey communities and how the effects of predator cannibalism on prey community composition are eventually reflected in ecosystem function will be fruitful next steps.
While our knowledge of the ecological consequences of predator cannibalism has gradually accumulated (Crumrine, 2010a;Persson et al., 2003;Rudolf, 2006Rudolf, , 2012Takatsu et al., 2017), community-level consequences of predator cannibalism are still unclear. Here, I showed that cannibalism among predators can reduce predation on small prey species but increase predation on large prey species simultaneously. Previous studies showed that the occurrence and strength of cannibalism among predators are highly dependent on surrounding environmental conditions, including the presence and absence of top predators, ambient temperature, and changes in water level (Crumrine, 2010b;Gillespie et al., 2020;Kishida, 2011;Sniegula et al., 2019). All this considered, cannibalism among predators could play a pivotal role in establishing links between biotic and abiotic factors and predator's effects on the prey community and resultant cascading effects on lower trophic levels. Investigating the context-dependency of cannibalism among predators further and how such context-dependent cannibalism among predators affects the nature and strength of top-down effects might be essential next steps in furthering our understanding of community-and ecosystem-level consequences of predator cannibalism.

F I G U R E 4
The numbers of (a) Ephemeroptera, (b) Diptera, and (c) Hemiptera remaining in the enclosures over the experimental period. Error bars denote standard error. The shaded area in the figures represents the period when cannibals, non-cannibals, frog tadpoles, and aquatic insects coexisted within enclosures. The total number of (d) emerging Ephemeroptera and (e) emerging Diptera was collected during the experiment. The thick horizontal bars represent the median, the box contains 50% of the data, and the whiskers indicate the range. Treatments not sharing the same lowercase letter were significantly different from each other (Tukey's HSD: p < .05). See Figure 2

DATA AVA I L A B I L I T Y S TAT E M E N T
Data available from the Dryad Digital Repository https://doi. org/10.5061/dryad.8cz8w 9gt5 (Takatsu 2022).